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Use  of  Continued  Fractions  in  High  Speed  Computing'' 

D„  Teiehroew 

University  of  North  Carolina 
and 

National  Bureau  of  Standards 

10  Introduction o In  the  course  of  carrying  out  computations  required 
for  numerical  solution  of  problems  it  is  frequently  necessary  to  have  avail- 
able the  value  of  one  or  more  functions  for  various  values  of  their  argu- 
ments o If  a high  speed  machine  is  being  used  it  is  usually  not  efficient 
to  look  up  tables  of  the  functions  outside  the  machine  or  store  them  in 
the  internal  memory  and  it  is  therefore  necessary  to  calculate  the  values 
of  the  function  -whenever  they  are  required e This  is  generally  done  either 
by  using  a rational  approximation  to  the  function  or  a finite  number  of 
terms  of  an  infinite  process G Usually  the  Infinite  process  used  Is  that 
of  the  power  series 0 The  purpose  of  this  paper3  is  to  show  that  another 
infinite  process,  a continued  fraction  expansion,  may,  in  some  cases,  be 
more  efficient e 

The  choice  of  which  infinite  process  to  use  may  depend  ons 
10  Properties  of  the  method  of  computing,  such  as 

(a)  the  number  of  orders  required  to  program  the  calculation 

(b)  rounding  error  involved 

(c)  the  magnitude  of  the  numbers  entering  into  the  calculations 

^This  work  was  performed  on  a National  Bureau  of  Standards  contract 
with  the  University  of  North  Carolina » 


. 

Sal' 

' 

. 


* 


2„  Properties  of  the  series , such  as, 

(d)  speed  of  convergence 

(e)  region  of  convergence 


3*  Properties  of  the  machine  being  used,  such  as 

(f)  amount  of  high  speed  storage  available 

(g)  what  factors  determine  the  total  time  spent  on  a problem «, 

Since  methods  of  computing  continued  fractions  are  not  as  well  known 

as  they  deserve  to  be,  the  following  section  will  be  devoted  to  examples  of 
three  methods  of  computation 0 In  section  3 the  speed  of  convergence  of 

some  particular  continued  fractions  will  be  compared  with  that  of  the  cor- 
responding power  series  and  in  section  h we  shall  discuss  some  situations 
in  which  continued  fractions  might  be  usedc 

20  Computing  Continued  Fractions 0 A continued  fraction  expansion  of 
the  function  f(x)  takes  the  form 


where  the  aDs  and  b^s  may  be  functions  of  xc  For  convenience  the  expansion 
is  usually  written  in  the  form 


(i) 


f(x)  = t»0  + aj_ 


'2 


(2) 


f(x)  - bo  + &2 


♦ 


The  following  are  examples  of  such  expansions 
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(3) 


x _ 
e = 1 


«ip 

Jii.  iV 

1 + 1=3+  ^“>1+  8»0 


(U)  In  x a x-1  X“(x°l)  X“(x-X)  2^(x-l)  2 (x-l)  3^(x~l) 


T“  + — T-  + 


+ p 


(?) 


arc  tan  x - x (x)^  (2x)“  (3x)^ 


1 + 


+ 7 + 


In  practice  a finite  number  of  terms  are  used  to  approximate  the  function^ 
the  so-called  "nth  approximant"  is  given  by 


(6)  fn(x)  = bo  + ^ Sg  an 

b-.  + b_2  * 000  bn 

Three  different  ways  of  computing  f (x)  are  given  by  the  following  methods 
I,  II  and  III 

Method  lo  The  obvious  way  to  compute  f (x) 3 if  n is  given5  is  to 
cany  out  the  successive  additions  and  divisions  indicated  by  the  form  of 
the  expansion  when  it  is  written  as  in  (l)  0 More  formally , the  method 
consists  of  calculating  in  sequence 


1 . « b . + 

n-i  n-i 


'n-i+1 


ml 


^ 0 


© 


n-a 


n-1  d. 


'n-i 


for  i - 0,  1 «««  n - le  Then  f (x)  - b + e,  • This  method  is  illus- 

33  n o X 

trated  in  Table  1 by  the  computation  of  arc  tan  1 using  10  terms  of 
expansion  (?)e 


i «e  f5  .1 

ea 
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Method  II,  Successive  approxiraants , for  consecutive  n,  may  be  caleu- 
lated  by  writing 

A 


f (x)  _ 

nv  7 B 


n 


n 


where  A s b A,  B b b.  + 

o o 1 o 1 


A _ = b _A  + a _A 
ml  ml  n ml  n- 


B - 1 B _ « b. 

o 11 


B s fe  _ B + a - B - 
ml  ml  n ml  n-1 


An  example  using  this  method  is  given  in  Table  2 where  expansion  (I4.)  is 
used  to  calculate  In  2 *3026 , 

Method  III,  The  nth  approximant  may  also  be  expressed  as  a sum,  i,e. 


n 

f (x)  58  b + 2 f ’ /I  f. 

nx  ' o ^ 2 1 


where 


ri  = b.  ,b.  : ' A " ,) 

1-1  1 1^  ^l-l' 


A * 1 * ^2 


1 


Table  3 gives  the  computation  of  In  2 O3026  by  this  method 
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Table  X Computation  of  arc  tan  1 using  ten  terras 
of  the  continued  fraction. 


b . 
n-i 

d . 

rt-i 

a 0 

n-x 

c . 

rt-i 

19 

19, 

81 

U.26315789 

17 

21 .26315789 

6U 

3.00990099 

15 

18 .00990099 

h.9 

2.72072567 

13 

15.7207256? 

36 

2.28997063 

11 

13 .28997063 

2? 

1.88111778 

9 

10.88111778 

16 

1.1*701*3717 

7 

8.1*701*3717 

9 

1.06251895 

5 

6.06251895 

k 

.65979175 

3 

3 .65979175 

1 

.27323959 

1 

1.27323959 

1 

.70539811* 

' 


- . 
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n 

0 

1 

2 

3 

h 

5 

6 

7 

8 

9 

10 

11 

12 

13 

1 h 

15 
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Table  2*  Computation  of  In  2 *3026  by  method  II, 


b 

n 

an 

A 

n 

Bn 

A /B 
n n 

0 

- 

0 

.1 

x 101 

- 

1 

1.3026 

.13026  X 101 

a 

x 101 

1.3026 

2 

1.3026 

.26052  x 101 

.33026 

x XO1 

.78883301*1 

3 

1.3026 

.951236676  x 101 

012101* 

2 

x 10 

.81*8530539 

h 

5.2101* 

.5162360112  x 102 

.6201*91*6? 

x 10' 

.831971*933 

5 

5.2101* 

.30768121*11*  x 103 

.3686580032 

x 103 

.831*598025 

6 

11.723U 

.21*51291571*  x Hr 

.293937871*1 

x 10^ 

.83391*8868 

7 

11.7231* 

5 

.2076611128  x Hr 

.21*8975761*2 

x IQ3 

0&3bp6L$5k 

8 

20.81*16 

.2172177286  x 10^ 

.2601*1*19673 

x XO6 

.831*03501*8 

9 

20.81*16 

.238775851*3  x 107 

.2862883031* 

x XO7 

.831*039852 

10 

32.565b 

.3095128077  x IQ8 

.3711012300 

x XO6 

.831*038751* 

11 

32.565b 

.1*182211*1*53  x 109 

.5011*lill389 

x XO9 

.831*038959 

12 

1*6.8936 

.61*70071*325  x 1010 

.7757520931 

■m10 
x 10 

.831*038913 

13 

1*6.8936 

.1037228751*  x XO12 

.121*3621523 

x 1012 

.831*038922 

111 

63.82?!* 

.1865088276  x 1013 

.2236212522 

x XO13 

.831*038919 

15 

63 ,8271* 

.31*59668557  x 10l!* 

.1*11*8090065 

x 10xit 

.831*038920 

- 


I 


■ 

' 


- 


' 
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3 
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10 

n 

12 

13 

iU 

15 
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Table  3,  Computation  of  In  2,3026  by  method  III 


rn 

1 + f 

n 

1.3026 

1.0 

.63.3 

.6055831*80 

.2171 

.883801*325 

.1*31*2 

.722675623 

.26032 

.81*1558660 

.39078 

.752522289 

.279128571 

.8261*21810 

.372171Ii29 

.761*77921+6 

.289U66667 

.81871*7283 

.361833333 

.7711*56089 

.29601*5*455 

.811*076311 

.355251*51*5 

.77567221*5 

.300600000 

.810920125 

.350700000 

.778579901* 

.30391+0000 

.80861*171*3 

f 

n 

1.3026 

1.3026 

-.513766959 

.78883301*1 

.0596971*99 

.81*853051*0 

-.016555608 

.831971*932 

.002623093 

.831*598025 

-.00061*9157 

.83391*8868 

.000112686 

.831*061551* 

- .000026506 

.831*03561*8 

.000001*801* 

.831*039852 

-.000001098 

.831*038751* 

.000000201* 

.831*038958 

-.00000001*6 

.831*038912 

.000000009 

.831*038921 

-.000000002 

.831*038919 

.000000000 

.831*038919 
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3®  Truncation  Error*  These  errors  arise  because  only  a finite  number 
of  terns  of  an  infinite  process  can  be  used*  For  a discussion  of  these 
errors  it  is  convenient  to  separate  functions  into  two  classes  on  the  basis 
of  the  order  of  magnitude  of  the  nth  term  in  their  power  series  expansions* 
If  the  order  of  magnitude  of  the  nth  term  is  xn/nl  the  function  belongs  to 

r>  . 

the  first  class;  if  it  is  x /n  the  function  belongs  to  the  second  class* 

The  reason  for  this  classification  is  that  functions  in  the  first  class  can 
generally  be  computed  adequately  from  their  power  series } while  for  functions 
of  the  second  class  an  unreasonable  number  of  terms  may  be  required. 
Heising^  states  that  in  computing  arc  tan  x by  means  of  a general  purpose 

board  on  the  6OI4.  Electronic  Calculating  Punch  the  maximum  error  after  990 

-7 

terms  are  used  is  still  11  x 10  . Fortunately  the  continued  fractions  of 

some  of  the  functions  in  the  second  class  converge  rapidly  enough  to  make 
them  practical  for  machine  computation® 

x 

Examples  of  functions  belonging  to  the  first  class  are  e , sin  x, 
cos  x5  sinh  x and  cosh  x and  examples  of  functions  belonging  to  the  second 
class  are  In  x>  arc  sin  x,  arc  sinh  x arc  tan  x and  arc  tanh  x*  To  com- 
pare the  convergence  of  the  power  series  and  the  continued  fractions  ex- 
pansions  of  functions  in  these  two  classes  we  have  chosen  e to  represent 
the  first  class  and  In  x and  arc  tan  x to  represent  the  second.  The  con- 
tinued fractions  expansions  are  given  by  (3)  3 (U)  and  (?)  and  the  power 


series  are 


' ■ . 

- ■ . • 

( 

. 

' 


• • • 


■'  . 


. . • - ■ c t.  . 


■■  ~ '1  . X r *U  ' • 


. • 


' OS 

■ 

- 

' • ■ ■ 
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(7) 


e 


1 + 


© © e 


(8) 

(9) 


In  x = 


(x  - 1)  (x-1)  ■ x (x—l) ‘ 

5 + —T— 


arc  tan  x - x 

T 


o © • 


Tables  U,  5 and  6 give  the  comparisons,  The  calculations  were  carried  out 
to  two  more  decimal  digits  than  were  required^  this  should  make  the  rounding 
error  negligible  (except  perhaps  for  the  computation  of  In  x for  x near 
zero)  , 


Table  U,  Truncation  error  in  e after  10  terms 

/ 


X 

Power  Series 

Continued  Fraction 

1 

2.7  x lO-8 

*67  x 10  ^ 

2 

61  x 10~6 

22  x 10“6 

3 

59  x 10-1* 

Uh  x 10  ^ 

Table  h gives  the  difference  between  e and  the  10th  approxiraant 
from  the  power  seri.es  and  from  the  continued  fraction  (in  the  power  series 
the  last  term  is  x"^/(L0?. . This  table  indicates  that  the  continued  fraction 
converges  a little  faster  than  the  power  series  but  the  difference  is  too 
small  to  be  of  practical  importance. 


‘ 
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Table  5o  Number  of  terns  required  to  compute  In  x t©  9 decimals 


X 

Power  Series 

Continued  Fraction 

<,0001 

95,000* 

550 

.0010 

11,500* 

315 

.0101 

i,35b* 

105 

,io5U 

lii3 

3U 

.5108 

25 

16 

.6931 

16 

11 

.9163 

8 

7 

2.3026 

=* 

16 

U.6052 

— 

2k 

6.9078 

30 

Approximate  values 

Table  5 gives  a comparison  ©f  the  number  ©f  terms  required  for  fixed 
accuracy®  For  the  power  series  the  number  of  terms ,<  n,  is  determined  by 
the  condition  that  the  absolute  value  of  the  (n  + l)th  term  in  (8)  shall 
be  less  than  9 <3132  x 10  The  continued  fractions  were  computed  by 

Method  III  and  n was  obtained  as  the  first  Integer  for  which 
(f^  0 0 0 ^n)  9 <»3132  x 10  10o  The  table  shows  that  not  only  does  the 

continued  fraction  expansion  converge  appreciably  faster  than  the  power 
series  but  it  also  converges  for  all  x>  Q,  while  the  power  series  con- 
verges only  for  0 < x ^ 2C  The  continued  fraction  could  be  used  t©  com- 
pute,  In  a reasonable  number  of  terms,  In  x for  01  <x  <10o 


11  - 


Table  6e  Number  of  terms  required  to  compute  arc  tan  x to  6 decimals 


x 

a 


o2 


oU 

o 6 
o? 
08 
3 
leO 
2 0Q 


Power  Series 
3 
U 

5 

6 
8 

11 

x5 


Continued  Fraction 


a 

5 

5 

7 

7 

7 

8 
8 

15 


Table  6 gives  the  corresponding  comparison  for  the  arc  tangent 0 For 

the  power  series , n,  is  determined  by  the  condition  that  the  absolute  value 

of  the  (n  + l)th  term  in  (9)  shall  be  less  than  10  while  for  the  continued 

fraction  n is  the  number  of  terms  required  to  reduce  the  truncation  error 
=6 

to  less  than  10  "0  The  continued  fraction  converges  faster  than  the  power 
series  and  also  converges  in  the  region  where  the  power  series  does  nots 
i0e0>  for  x greater  than  lc 


. 


' 


f .1  f-  ft  ^ i | '{i||  m ii  9 j 
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Ue  Machine  Computation 0 Computing  machines  differ  greatly  in  many 
aspects  o In  deciding  what  method  should  be  used  to  compute  values  of 
functions  the  amount  of  high  speed  storage  available  is  one  of  the  most  im- 
portant factors  o 

Machines  which  have  enough  high  speed  storage  to  accommodate  arbitrary 
coefficients  can  probably  compute  values  of  many  functions  most  rapidly  by 
using  polynomial  approximations  with  a fixed  number  of  terms  0 If,  however, 
the  majority  of  the  arguments  fall  in  the  region  in  which  only  a few  terms 
of  an  infinite  process  are  required,  or  if  varying  degrees  of  accuracy  are 
required  from  one  computation  to  another,  the  infinite  process  may  be  more 
efficient.  In  general  the  computation  of  a term  in  a power  series  reqiires 
fewer  operations  than  the  computation  of  a term  in  the  continued  fraction 0 
However  for  functions  such  as  In  x and  arc  tan  x more  terns  of  the  power 
series  are  required.  In  either  case  it  is  desirable  to  compute  only  as 
many  terms  as  are  required  to  achieve  the  required  accuracy.  Therefore, 
for  computing  continued  fractions,  method  I would  not  be  suitable,  but 
method  II  or  III  could  be  used.  It  is  immediately  evident  from  Table  2 
that  method  II  cannot,  in  general,  be  used  without  a floating  decimal 

point*  because  A and  B increase  too  rapidly 0 Furthermore,  since  method  II 

* n n 

requires  as  many  operations  and  more  storage  than  method  III  the  latter  is 
more  suitable  for  machine  calculations. 

Table  7 gives  the  calculations  required  In  computing  successive  ap= 
proximants  to  the  continued  fraction,  the  values  required  from  storage 
and  the  initial  values,  if  method  II  is  used.  If  a and  are  simple 
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functions  of  n and  x which  can  be  computed  as  required  the  first  calculation 
requires  the  storage  of  only  n and  xc  For  example y for  expansion  (k) 


T 


n 


Table  ? 


Calculation 


rn  VVlbn 


Required  from  storage 

a , b , , b 
n’  in-l’  n 


Initial  values 


20  1 + 


n 


nrl 


Pt 


2 / n~l 


’ n n-1  t 1J 


x 


where  [y]  denotes  the  greatest  integer  ^ y0  Even  with  this  simplification 9 
five  quantities  must  be  stored e 

This  fact  makes  this  method  of  computation  impossible  for  machines 
with  very  small  internal  memories3  such  as  the  IBM  60li  Electronic  Calculate 
ing  Punciie  The  computation  of  power  series  requires  very  little  storage 
and  can  be  carried  out  on  such  machines  0 However  for  functions  of  the 
second  class  where  too  many  terms  of  the  power  series  are  required 3 con- 
tinued fractions  can  be  used  if  the  maximum  number  of  terms  of  the  eon- 
tinned  fraction  expansion  that  will  be  needed  to  give  the  desired  accuracy 
is  determined  in  advance  and  method  I is  used  to  compute  the  value  of  the 


. 


, 


- 1U  - 


expansion e If  the  coefficients  can  be  generated  in  a simple  manner  only 

1 

three  numbers  have  to  be  stored  at  any  one  t±me0 

So  Bibliography e The  general  theory  of  continued  fractions  is  given 
1 2 

by  Perron  and  Wall  „ Both  of  these  books  have  extensive  bibliographies,, 

The  use  of  continued  fractions  in  interpolation  is  discussed  in 
3 oo  li  5 & 

Milne-Thomp son  , Noriund  and  Lane  0 Muir  gives  some  general  methods  for 

oo  °7 

transforming  infinite  series  into  continued  fractions  „ Muller  and 

8 9 

Aroian  and  Burgess  give  some  examples  of  the  use  of  continued  fractions 
in  evaluating  integrals 0 


The  use  of  continued  fractions  in  computing  In  x,  arc  tan  x and 
arc  sin  x on  the  Model  IX  Card  Programmed  Calculator  will  be  discussed 
in  a forthcoming  Bureau  of  Standards  report 0 
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